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Abstract. An approach to modeling fracture incorporating interfacial mechanics is applied to the example 
of a curvilinear plane strain crack. The classical Neumann boundary condition is augmented with curvature- 
dependent surface tension. It is shown that the considered model eliminates the integrable crack-tip stress and 
strain singularities of order 1/2 present in the classical linear fracture mechanics solutions, and also leads to the 
sharp crack opening that is consistent with empirical observations. Unlike for the case of a straight crack, for 
a general curvilinear crack some components of the stresses and the derivatives of the displacements may still 
possess weaker singularities of a logarithmic type. Generalizations of the present study that lead to complete 
removal of all crack-tip singularities, including logarithmic, are the subject of a future paper. 

Key words. Fracture mechanics, curvilinear crack, surface tension, complex potentials, singular intcgro- 
diffcrential equations. 

AMS subject classifications. 74B20, 74G70, 74R10. 

1. Introduction. A better understanding of fracture processes is of utmost importance 
for applications. The standard approach to study the behavior and propagation of fracture is 
based within the framework of linear elastic fracture mechanics (LEFM). This approach has been 
successfully applied to countless problems and has a vast literature. However, LEFM contains 
a well-known internal inconsistency. While it is based on the assumption that the stresses and 
the strains remain small everywhere in the body, LEFM predicts that the stresses and the 
strains possess an integrable power singularity of the order 1/2 at the crack tips. Moreover, in 
the case of an interface crack between two dissimilar materials there is an additional oscillating 
singularity which results in the non-physical interpenetration of the crack faces near the crack 
tips. 

Several models have been proposed to eliminate this internal inconsistency. One of the 
most common approaches, studied by many authors, is to introduce two-dimensional cohesive 
zones or three-dimensional process zones near the crack tips. Despite the obvious advantages 
of this approach it has shortcomings of both theoretical and practical nature, such as difficulty 
in specifying physically valid constitutive response functions. 

Since the fracture occurs within nano- or molecular scale processes, it has been argued that 
it is impossible to describe the fracture effectively within the context of continuum mechanics. 
Thus, there is a growing literature dedicated to the modeling of fracture using atomistic and 
lattice based approaches [H [51 |3J [SJ [TI>J HH [23 US]- The accuracy of these methods largely 
depends on the precise description of intermolecular forces which is difficult to do for liquids 
and solids [15] . This approach may also present some computational challenges. Several atom- 
to-continuum models have been proposed as well. One of the most extensively studied methods 
of this type, in the context of finite element method (FEM) approximations to continuum 
models, is the quasi-continuum method introduced by Tadmor et al in 1996 [3D]- Based on an 
atomistic view of material behavior, its continuum aspect comes from the fact that the FEM 
is based on energy minimization. A different type of atom-to-continuum modeling is a recently 
proposed approach by Xiao and Belytschko in [31] which involves the introduction of bridging 
domains between regions modeled using bulk (continuum) descriptions of material behavior and 
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regions modeled using atomistic descriptions of material behavior. Both of these approaches 
involve adjustable parameters, such as choosing the size and location of domains over which 
potentials acting at different length scales are blended, that one needs to ht for every particular 
application. 

In a series of works [Tl IIQ1 [III H21 H3 [2B] the modeling of fracture is based on the introduction 
of the concept of configurational forces. The configurational forces were first introduced as 
derivatives of energy by Eshelby [5] and were further studied by Gurtin [5] . The configurational 
forces are described in j9] as the intrinsic properties of a body's material micro-structure which 
have their own force balance. 

Another continuum approach, peridynamics [T31 HH addresses the modeling of the 
spontaneous formation and dynamic propagation of a fracture, for which the classical continuum 
mechanics approach is ill-suited as well. In peridynamics, material points are assumed to be 
connected via elastic (linear or nonlinear) bonds which have a critical relative elongation. This 
model does not distinguish between the points in a body where a discontinuity in displacement 
or its derivatives may be located. The essence of the mathematical aspect of this nonlocal 
model is that integration rather than differentiation is used to compute the force on a material 
particle. 

One of the most convenient and successful modifications to LEFM incorporates surface 
excess properties on the crack boundary into the modeling the fracture. It is well-known that 
the molecules on the boundary between two different phases are subjected to different forces 
than molecules in the bulk of the material. This model accounts for this fact by introducing 
a dividing surface endowed with surface excess properties |25) . This theory has been first 
applied to fracture mechanics in |20j and developed in more details on the example of a straight 
mode-I crack in [35]. The advantages of this model are that, first, it relies on the physically 
valid assumption that the forces which act on the boundary of the object are different from 
the forces in the bulk and a curvature-dependent surface tension is present along the interface 
between two phases. It has been observed also that the behavior of materials is in accordance 
with linear elasticity everywhere except very small neighborhoods around the crack tips which 
provides additional justification for using the equations of linear elasticity to describe the bulk 
of material. Hence, since the standard equations of linear elasticity are assumed for the bulk 
of the material and only the boundary conditions arc changed, it follows that a wide variety 
of mathematical approaches and techniques used in linear elasticity are still applicable to this 
model. Finally, this model predicts bounded stresses and strains in the case of a straight 
mode-I crack, thus removing the inconsistency present in LEFM. This model also bypasses the 
introduction of artificial conditions on the boundary and avoids ad hoc choices of parameters 
that are present in most other models. 

While the body of research discussed above is dedicated to modification of the LEFM in 
order to remove the power singularity of the order 1/2 at the crack tips, it is necessary to note 
that most of these attempts, with rare exception, consider straight cracks only. For practical 
applications, however, it is important to study curvilinear cracks. In this paper we extend the 
curvature dependent surface tension model [22) to the curvilinear smooth plane strain crack of 
essentially arbitrary shape. The most significant conclusion of the paper is that the analysis 
of a curvilinear (non-zero curvature) crack is significantly different from the case of a straight 
crack. Even though just as in the case of the straight crack, the stresses and strains do not 
have an integrable power singularity of the order 1/2, thus removing the main inconsistency 
of LEFM, the shear stress and one of the derivatives of displacements may possess a weaker 
logarithmic singularity while the tensile stress and the derivative of the other displacement 
remain bounded. Thus, the curvature of the crack plays an important role in the singular 
behavior of the stresses and strains. This suggests that controlling mean curvature alone may 



Modeling of a curvilinear planar crack. 



3 



not be enough to force the stresses and strains to be bounded. We conjecture that the present 
theory must be extended through the addition of higher surface gradient effects in the model 
for stress-deformation behavior within the fracture surfaces in order to remove the logarithmic 
singularities predicted by the curvature dependent surface tension model considered here. Such 
considerations are beyond the scope of the present contribution and are the subject of an ongoing 
research effort. 

The complex analysis approach is used for the solution of the stated problem. Functions of 
complex variable have been first applied to the solution of two-dimensional problems of elasticity 
by Kolosov in 1909. The methods of complex analysis have been further developed by Muskhe- 
lishvili |19) who inspired a large school of research in this area. The complex analysis approach 
is an extremely powerful technique and can be applied to a large number of problems of two- 
dimensional elasticity. The main idea consists in expressing the stresses and the displacements 
in elastic body through two analytic functions, thus, reducing the problem of elasticity from 
the solution of partial differential equations to the solution of boundary problem for these two 
analytic functions. This boundary problem is often reduced further either to Riemann-Hilbert 
type of problem or to integral equations of some kind. Some excellent applications of methods 
of complex analysis in elasticity can be found in [4], [19], [27] and others. 

The solution of the problem proceeds in the following way. First, the boundary condition 
of the type [35] is stated. The curvature of the crack dictates a slight modification to the 
form of the surface tension considered in [22]. In particular, it is reasonable to assume zero 
surface tension in the unloaded configuration of the body. Next, this condition is linearized 
and written in terms of the tensile and the shear stresses and derivatives of the displacements. 
Muskhelishvili's formulas [H] are applied then to reduce this condition to the boundary value 
problem for two analytic functions (complex potentials) and their derivatives. On the next 
step, the boundary value problem is reduced to the system of two Cauchy-type singular integro- 
differential equations with the help of Savruk's integral representations of complex potentials 
[21] . This system is further reduced to the system of two weakly singular integral equations 
in the spirit of [18]. It is necessary to note that the mathematical techniques used to solve 
the problem are significantly different from those in |22) . Finally, the numerical procedure is 
developed to solve the system of singular integro-differential equations and the numerical results 
are presented. 

2. Model with a curvature dependent surface tension. Consider the differential 
and jump momentum conditions in the deformed configuration in the absence of inertial and 
gravitational effects [23] : 

div(T) = 0, (2.1) 
grad (c) 7 + 27i/n+[T]n = 0, (2.2) 

where T is the Cauchy stress tensor, n is the unit normal to the fracture surface £ pointing into 
the bulk material, H = — ^div^n is the mean curvature, subscript . . .iq is used to emphasize 
that the quantity is associated with a surface of the crack, grad^ and div^) denote the surface 
gradient and the surface divergence correspondingly with respect to the given parametrization 
of the surface and can be derived using the formulas in [24] , and the double brackets [. . .] denote 
the jump of the quantity enclosed across the boundary of the crack. 

Assume that 7, which will be called surface tension, depends linearly on the difference 
between the curvature of the deformed crack and the curvature of the crack in the unloaded 
configuration: 



7 = 7i(div (c) n - div (Co) n ), 



(2.3) 
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where the subindex "0" denotes the parameters in the initial configuration of the crack. 




Fig. 2.1. An infinite plate S with a crack L. 



Since the crack L in an undeformed configuration is a Jordan arc (fig. 12. ip . it is natural to 
take the parametric equations of the crack L to be t(s) = Xi(s) + ix2(s), where the parameter 
s is an arc length, s E [0,/]. Assume that the function t(s) has continuous derivatives up to 
the fourth order. This assumption is made for simplicity and can be somewhat relaxed. Since 
the boundary of the crack is free from the external stresses and, hence, the jump of the stresses 
on each of the boundaries of the crack is simply equal to the stresses in the cracked plate, it is 
possible to rewrite the condition (|2.2p in terms of the tensile and shear stress in the following 
form: 



c„ + ir n = i(grad (c) 7 - 7ndiv( C) n)f'(s), 



(2.4) 



where a n and r„ are the tensile and shear stresses on the crack boundary. Furthermore, it is 
easy to see that in the deformed configuration of the crack the quantity div^n can be expressed 
as 



div(£)n 



A2A1 — X\X2 



(2.5) 



where Xi(s) + 1X2(3) = (xi(s) +ui(s)) + i(x2{s) + U2(s)) are the equations of the crack in the 
deformed state and (u\ + iu2)(s) are the displacements. Linearizing the condition (|2.5j) under 
the assumption that the derivatives of the displacements are small, one obtains 



div (c) n = k (s) + (x 2 (s) - 3xi(s)k (s)) 



dui 



(xi(s) + 3± 2 (s)k (s)) 



du 2 
ds 



X2{s) 



d 2 ui 
~ds^ 



±i(s) 



d 2 u 2 
~ds r ' 



(2.6) 



where kq(s) is the curvature of the crack in the undeformed configuration. Substituting (|2.6I) 
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into (|2.4I) leads to the following boundary condition 

(a n + «r„) ± (s) = mi(s)^(«i + W2) ± (s) + m 2 {s)-^{u 1 - iu 2 ) ± (s) + m 3 (s)-^ (ui + ?U2) ± (s) 
+ "i 4 (s)^-(ui - w 2 ) ± (s) +«7ilm^- (t'(s)(ui +iu 2 ) ± {sfj , s G [0,1], 

(2.7) 

mi (s) = -O.571 (t"'(s) + 2it"(s)K (s) + 3i¥{s)K' {s) + 3^(s)/^(s)), 
m 2 (s) = O.571 (t"'(s) — iit" (s)ko(s) — 3z£'(s)/Cq(s) — 3i'(s)/c (s)), 



m 3 (s) = -2i7ii'(s)Ko(s), m 4 (s) = -i7ii'(s)«o(s)- 

where "+" corresponds to the stresses and the displacements on the left-hand side of the crack 
with respect to the direction chosen, and "— " corresponds to the same functions on the right- 
hand side of the crack. Thus, the equation (|2.7[) is valid on the both surfaces of the crack. 

3. Complex potentials and their integral representations. It is well known [T§] 
that the tensile a n and the shear r„ components of the stress vector acting on the tangent 
line to the curve L from the side of the positive normal and the derivative of the displacement 
vector u± + iu 2 at the points of the curve L can be obtained from the Muskhelishvili's complex 
potentials $(z), &(z) by the following formulas |21j : 

(a n +ir n )(t) = ^ t )+W) + ^(tW)+Wj), (3-1) 

2/il(ui + iu 2 )(t) = «*(t) - ¥(t) - ^-(tWiT) + ¥(£)), i G L. (3.2) 
or at 

Here &(z) and ^(z) are analytic functions outside of the crack, /i is the shear modulus of the 
material and k — (3 — + v) for the case of plane stress and n = 3 — 4is for plane strain, 

where v is the Poisson's ratio. 

Consider the infinite plate S with a crack L under the action of the given at infinity stresses. 
Assume that the stresses o„ + ir n and the derivatives of the displacements 4z(ui + iu 2 ) in the 
plate have jump discontinuities on the curve L. Then the complex potentials <&(z) and ^S(z) 
can be taken in the form I2T1: 




2tt Jl\ t-z (t- zf 

zeS, r = K° + of )/4, T' = (o 2 °° - o-)e- 2lQ /2, 

where cr^* and o|° are the principal tensile stresses at infinity acting in the directions which 
constitute the angles a and a + it/ 2 with the positive direction of cc-axis correspondingly, and 
g'(t) and q(t) are the unknown functions which can be expressed through the jump of the dis- 
placements u\ + iu 2 and the jump of the stresses a n + ir n on the line L by the formulas: 

g'(t) = ((mi + iu 2 )+{t) - ( Ul + iu 2 )-(t)) , t G L, (3.4) 

q(t) = ((o„ + ir n )+(t) - (o„ + »t„)-(*)) /2, * G L. 
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Hence, the stressed state of the cracked plate is described by the complex potentials 
which contain two unknown functions g'(t), q(i) defined on the crack L. We will look for these 
functions in the class of functions satisfying the Holder condition on the curve L. This choice 
guarantees [8] the existence of all principal and limit values of the integrals of the Cauchy type 
in the formulas (|3 .3|) . Also assume that the formula (I3.2|) can be differentiated by s twice. This 
will be proved rigorously later. 

Using (|3.1|) . (|3.2[) and (|3.3|) it is possible to express the stresses (er n + iT„) and the derivatives 
of the displacements -r-(u\ + iu%) in the form: 



(<r„ + ir„) ± (s ) = ±q(s) + 



2tt 



\-kx(s,s ) )g'(s)ds + / k 2 (s, s )g'(s)ds 



-2i 



— h k 3 (s, s ) ) q(s)ds + 2i / k 2 (s, s Q )q(s)ds 

o \ s-«o / Jo 



2Rer + r'(t'(s )) 2 ; (3.5) 



d 



2/x— (ui + m 2 ) ± (so) = t'(s )O ± (s ) + (kT - T)t'(so) - r't'(s ), 



ds 

where the functions fi^So) have the following form 



(3.6) 



" (*o) = ± 2 V («o) + 2?r ; 



« - 1 

o V s — s o 



+ fe 4 (s, s ) ) g'(s)ds 



(3.7) 



fc 2 (s, so) g'(s)ds - 2i 



2 k 



+ Kk±(s, sq) ) q(s)ds — 2i / fc 2 (s, so)g(s)ds 



where gr'(s) = g'(t(s)), q(s) = q(t(s)) and kj(s,so), j = 1, 2, 3, 4, denote the regular kernels: 

f(«) , 



fci(s, s ) 
k>(s, s ) 
h(s,s ) 



+ 



s-s ' t{s)-t(so) t(s)-t(so)t'(so)' 
F(s) {t(s)-t(s ))t'(s)t'(so) 



ki(s, s ) = - 



t(s)-t(s ) (t(s)-t(so)) 2 t'isoY 

K-l | t'(s) Kt'(s) t'(s Q ) 

s-s t(s)-t(s ) t(s) - t(s ) *'( s o) ' 

K-l Kt'{s) t'(s) t'(s ) 



(3.8) 



+ 



s - s t{s) - i(s ) t(s) - t{s ) t'{s ) ' 
It is possible to rewrite the boundary conditions (|2.7p in the following form: 



(<7„ +ir„) ± (s ) = -r^-Ko(so) 



2/i 

ill — 
2[i dso 



Ko(s )Refi ± (so) - Im— fi ± (s ) 
ds 



K (so)Refi ± (s ) - Im— fi ± (s ) 
dso 



/(so), 



(3.9) 



where the functions f2 (so) have been introduced above, and 



/(so) - li^Y 1 mi(« )((KT - r)t'(so) - r't'(so)) + m 2 (s )(( K r - T)t'(so) - T't'(s )) 



+ m 3 (s )(( K r - T)t"(s Q ) - rt"(so)) + m 4 (s )(( K r - T)t"(s ) - T't"(s Q )) 

+ 2i im{t'(a )((Kr - r)*"'(s ) - r'i'"(s ))}l -2R e r-r' {t'{s )) 2 . 
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The boundary conditions (|3.9p must be combined with an additional condition providing 
that the displacements are single- valued on the crack L which can be expressed in the simple 
form as: 



g'{s)t'{s)ds = 0. 



(3.10) 



Observe that since (a' rl -|-iT„) ± (so) and fi^So) are allowed to have at most integrable power 
singularities at the crack tips, it follows immediately that Ko(so) Re f2 ± (so) — im af^^ ± ( s o) must 
be bounded at the the crack tips, and hence, the stress <r„ and the derivative of the displacement 
duf/ds are bounded. Additionally, boundedness of Ko(so) Re r2 ± (so) — Im^fJ^So) makes it 
possible [18] to differentiate the Cauchy singular integrals in (|3.9[) . Substitute Q3.5p . (|3.7[) into 
(|3.9j) to obtain the following system of equations: 



1 

2^ 



(g'(s) + g'(s))ds (k-1) f l (q(s) - q(s))ds 



s - s 



2ni 



s - s 



7i 



+ J^^oiso) < «o(so) 



2tt 



k - 1 f l (g'(s)+g'(s))ds + 2k f l (g(s) - q(s))ds 
o s - s iri Jo 



s - s 



ds 



2vr 



s - so 



k-1 [ l (g'(s)-g'(s))ds + 2k [ l (q(s) + q(a))da 







s - so 



(3.11) 



ReMx(g',q)(s ), s € [0,/]; 



2iri Jo 



(g'(s)-g'(s))ds , f (q(s) + q(s))d, 



s - s 



. d 



2tt 



2tt 



s - so 



k-1 f l (g' (s)+g'(s))ds + 2k f l (?(«) ~ l( s )) ds 
s — so Tti 



(i 



s - s 



ds 



K-1 f l (g'(a)-g'(s))ds 2k f l (q(s) + q{s))ds 



2vr 



s - s 



s - s 



(3.12) 



= ImM 1 (j',9)(s ), 8 € [0,1]; 
where Mi(g', q)(so) is a Fredholm integral operator of the first kind 
1 



M 1 (g',q)(s ) = 



2tt Jo 



ki(s, so)g'(s) + k 2 (s, s Q )g'(s) - 2ik 3 (s, s )q{s) + 2ik 2 (s, s )q(s) ds 



^( K oOo) +in' {so))Re 



1 



1 

2Wo 



(fc 4 (s, so)g'(s) - k 2 (s, s )g'(s))ds 



H : / (Kki(s,s )q(s) + k 2 (s,s )q(s))di 



(3.13) 



7i 



K o(so) / (ik 4 .2{s, so)g'(s) - ik 2}2 (s, s )g'(s) + 2nik 1 . 2 (s, so)q(s) + 2k 2 . 2 (s, s Q )q(s)^j ds 



i — 1m 



2/j, 



/ (k4, 22 (s,s )g'{s) - k 2 ^ 22 (s, s )g'(s))ds 



1 

iri 



(Kk 1:22 (s, so)q(s) + k 2>22 (s, s )q(s))ds 



+ (k+1)/(«o), 
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where kj^is, So) and fej- 22(^1 s o) denote the first and the second partial derivatives of the kernels 
(|3.8[) by the variable Sq. Here and further, assume for simplicity that Ko(sq) 7= everywhere on 
L. 

Similarly, by substituting (|3. 41) into (|3.9p we obtain the following relations between the 
functions g'{s) and g(s): 



q(s ) + q{s ) = - — k (s ) K (s )(g (s ) - s'(s )) + % (so) + 5"(«o)) 
7i d 



(3.14) 



g(s ) - 5(s ) = j^-j- Ko(s )(g'(s ) - g'(s )) + i{g"(s ) + g"{s )) 



ds 



(q(so) +9o(so)) 



Observe that the expression in the figure brackets {. . .} is the same in both of the equations 
(13. lip , (|3.12p . Solve the equation p. lip for this expression and substitute into the equation 
(|3.12p . Then, after differentiating and inverting the singular integrals, we arrive at the following 
regularized system of integro-differential equations: 



7i(«-l) 



2k (s ) 



8(i 



«o(s ) I (g'(s ) +g'(so)) + i 



( K- 1 



\2k (s ) 2/x 



4 M 



K — 1 
! 

2 
1 



(fl"(*o) - 5"(so)) + 2k(q'(s ) + g'(so)) 



/ 1 



^s (?-s )7ry V K o(s) 



ReM 1 (« 7 ', g )( s )+M 2 ( ff ', (? )( s ; 



i^G/(so)-?R) 

K oQo) 
2k (so) 



K - 1 



Ko(so)(?(so) + g(so)) 



(g'(s )+gi{s )) + 



i{n - 1) Kq(5q) 
2 k (s ) 



(«Oo) - g(s )) 



i(K-l) , 



(q (s ) - g'Oo)) 



1 



1 



Vso(Z - so) ^ Jo 



tt" K (so) (g(so) - g(so)) 



y s(l — s)ds 
s - s 



(3.15) 



K (s)lmM 1 (g',q)(s) - K (s)— , 

as \kq(s) 



ReM 1 ( 5 ',g)(s; 



y s(Z — s)ds 



Do 



s - s 



TTV^oC - s o) 



s G [0,1]; 



(3.16) 
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where the regular kernels M 2 (g' ,q)(s) and M 3 (g' ,q)(s) have the following form: 
M 2 (</,<?) (s ) 



1 



1 



1 



+ 



+ 



27T J \K (s) K (so) 
i(K- 1) 



(g'(s)+g'(s))ds 



1 



1 



2tt J \k (s) kq(so) 



s - s 
(q(s) - q(s))ds 
s - s 



2L 



K-l 

I 



(kI(s) - Kq(so)) 



(</(s)+ 5 '(s))ds 



— / ( K o(s) - «o(so)) 



(I 



s - s 

(q(s) - q(s))ds 
s - s 



M 3 (g',q)(so) = — I (k (s) - «o(« )) 



2vr 



(«o(s) - Ko(«o)) 



s - s 
(g(s) +g(s))ds 
s - s 



2tt 



«o(s) Ko(so)/ 



s - s 



«(«-!) /* / K o( s ) K o( s o) \ (<?(s) - g(s))ds 



2tt 



s - s 



/o V K o(s) Ko(so), 

The constants Z?i and D 2 in the equations (|3.15[) and p,16[) can be specified by integrating 
these equations and using the additional condition p.lOp : 



D, = 



1 



o I 2k (s ) 

K-l 



{g'{s Q )+g'(s Q )) + 



i(n - 1) 
2k (s ) 



(q(sq) -q(so)) 



-t-ko(so) 
4/x 



(9' (so) +g'(s )) ~ 2iK(q(s ) - q(so)) 



(3.17) 



2i 

4 /t 



i^-(9"(so) ~ 9"(s )) + 2K(q'(s ) + q'(so)) 



dso + Re g'(so)t'(so)ds ; 



Do 



-i^l(g'(s )-7W)) 



K-l 



Ko(s )(q(s ) + q(s )) 



k'o(sq) 
2k (s ) 



(9'(so)+g'(so)) 



i(k- 1) k'q(sq) 
2 kq(s ) 



(q(sq) - q(s )) 



(q'(so) ~ q'(so)) > ds + Im / g'(s )t'(s )ds . (3.18) 



- (9"(so)+g"(s )) 

Taking the constants D\ and D 2 in the form (|3.17j) allows to satisfy the conditions p.!0|) 
automatically. Note that the system p,15[) and (|3.16[) has a Holder continuous solution g'(s), 
q(s) only if the following conditions are satisfied at the ends of the crack: 



i^r±(9'(s)-g>(s)) + 2 K ( q (s) + q(s)) 



l(g'(s)+g'(s)) + l -^ J ^-(q(s)-q(s)) 



= 0; 



(3.19) 



s=CU 



= 0. 



=0./ 
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The physical meaning of the conditions (|3.19[) is that the tensile stress a n and the derivative of 
the displacement G?u 2 /(is are bounded at the tips of the crack. 

Finally, we will reduce the system of two regularized integro-differential equations (I3.15[) . 
(|3.16|) to the system of two weakly singular Fredholm integral equations. In the spirit of [18], 
denote the highest order derivatives in the equations (|3. 15|) . (j3. 16[) as new unknown functions: 



<p(s ) = i(g"(s ) - g"(s )); ip(s ) = q"{s ) + q"(s ). (3.20) 
Integrating the formulas (13.201) and making use of the formulas (I3.14j) obtain: 



where 



q'(so) + q'(so) = / uJo(s,so)i/j(s)ds + Ci; 



q(s )+q{s )= I 0Ji(s 1 s )^(s)ds + Cis + C 2 ; 
I 



i(9'i s o) - 9'(so)) = / uJo(s,s )<p(s)ds + C 3 ; 



V i 



■-2 



9 (so) +3" (so) = — — p-r / uj 1 {s,s )^(s)ds + C 1 s a + Q 

71 k (s ) yj j 

+ n {s ) / w (s, s )ip(s)ds + C 3 ; (3-21) 



4/i f w (s,s ) 



g'(so) + g'(so) = — / fr - / ui(si,s)ip(s 1 )ds 1 + C 1 s + C 2 ) ds 

7i Jo «o(s) \Jo J 

+ / cj (s, sq)ko(s) / w (si,s)^(si)dsi + C 3 ) ds + C 4 ; 



1 ,/ x , 2k (so) 



i(q'(s ) - q'(s )) = ^^(so) + / w o(s, s Ms)ds + C 



Ko(So) K-l{so 

+ d / ^so)\ ( f l Ui{StSo)i)is)ds + CiSQ + C2 



ds \ K o( s o) 

i(<l{ s o) ~ q( s o)) = t-^ ( / w (s, s )ip(s)ds + d 

K0{S0) \ Jo 



K o( s o) 



wi(s, s )ip(s)ds + CiSq + C-2 



Wo(s,so) = |q ifs^[o'so]' ' w i( s ' s o) = y w (s, si)w (si,s )dsi. 

The four constants Ci, C*2, C3 and C4 can be found from the four additional conditions 
(|3T9]1 . Substituting the formulas (|3~2"TT) into the equations ([3T5]) . (l3~16l) we obtain a system 
of two weakly singular Fredholm equations of the second kind with respect to the functions 
tp(s), ip{s). This system has a unique solution with the exception of the discrete spectrum of 
the Fredholm operators. It also follows that the functions f(s) and ip( s ) can have at most 
integrable power singularities of the order 1/2 at the crack tips. Hence, the functions g'(s) and 
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q(s) are bounded at those points. Thus, the stresses and the displacements can have at most 
logarithmic singularities at the crack tips. On the other hand, the conditions (|3.19p provide 
that the stress and the derivative of the displacement duf / ds are bounded at the end points 
of the crack. Thus, only the functions and duf /ds may possess logarithmic singularities at 
the crack tips. 

4. Straight Crack. While zero curvatures ko(so) = have been explicitly excluded from 
the consideration until now, it is possible to study a straight crack ko(so) = using the same 
methods as above. The equations (|3.11[) . (|3.12j) in this particular case become: 

1 f l (g'(s) + g 7 {sj)ds i(n - 1) f l (q(s) - qjs))ds 



2w Jq s — Sq 2tt Jq s — Sq 

= Re Mi (5', q)(so), s e[0,l}; (4.1) 



1 /' {g\s)-g'(s))ds 7 i(«-l) f l (g"'(*) - g"'( S ))ds , 
7T- 5— / =ImMi(g ,q){s ), so G [0,1 ; 

where the expression for the regular integral operator Mi(g' , q)(so) can be obtained from the 
formula (|3.13l) by substituting «o(so) = 0. Furthermore, instead of the equations f|3.14[) we 
obtain: 



q(s ) + q(s ) =0, s g [0,1]; (4.2) 



'■('y(-^u) -'/(-..)) : -^(5"'(so)+5"'(so)), s G[0,Z]. 



The conditions p.!9p become in this particular case: 



5 '(0)-<7'(0)=0, g"(0)-g"(0)=0; (4.3) 



g'(l)-g'(l) = 0, g"(l)-g"(l) = 0. 

Just like in the previous section, it is possible to invert the singular integrals in the equations 
(|4.ip to obtain the system of two regular integro-differential equations. Analogously, introduce 
two new unknown functions: 



<p{a ) = g'"(s ) + g"'(s ), 1>{a ) = i(g"'(s ) - g"'(s )), s G [0, l], (4.4) 

and integrate them similarly to (|3.2ip . Substituting (|4.4j) into the system of integro-differential 
equations (|4.ip we obtain a regularized system of two weakly singular Fredholm equations 
of the second kind. Similarly, it follows that g'(s) and q(s) are bounded at the crack tips. 
Furthermore, both the normal and the shear stresses cr^ and rj^ and the derivatives of the 
displacements duf /ds, d 2 uf /ds 2 are bounded at the crack tips, which completely corroborates 
the results of [55] which were obtained there using different techniques. 

5. Numerical solution of the system (|3.1ip . (|3. 12|) . The system of two weakly sin- 
gular Fredholm integral equations obtained by regularization of the system Q3.15p . (|3.16p in 
general has a very complicated form and is inconvenient for numerical solution. Hence, for 
practical purposes it is easier to solve the initial system of integro-differential equations p. lip . 
(|3.12p . It is well known that the systems of singular integral and integro-differential equations 
of this kind are notoriously difficult to solve numerically with a good accuracy and a reasonable 
computational time. From the authors' point of view, one of the best ways to solve (|3.1ip , (|3.12l) 
is to take a Taylor polynomial approximations of the unknown function g'(s) in the following 
form: 

N 

g'(s) = J2(gi + igl)(s-l/2) k , (5.1) 

k=0 
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Fig. 5.1. Graphs of Reg' and Imp' for N = 16, N = 20 and N = 30. 
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Fig. 5.2. Graphs of stresses o„ , t„ and derivatives of displacements du^/ds, du^ /ds for <tJ° = 1, cr^° = 0. 



where g\ and g\ are real coefficients. Using the formulas (|3.14[) . one obtains the following 
expressions for the real and imaginary parts of the function q(s): 



N 



q(s) + «(*) = ^-k (s) [kgiis - l/2) k ~ x + g 2 k K (s)(s - l/2) k ] , 

" k=0 



(5.2) 



i(q(s) q(s)) = -£ £ W fc " ~ + S*(*«o(«)(* - VZ)*" 1 + «{,(*)(« - Z/2) fc )] 



JV 



2//. 
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The following formula is used for numerical evaluation of the integrals in the equations p. 111) . 
(E32): 



f l I N 

/ K(S, S Q )(f)(s)ds = — — - ^ K(T k ,Sj)(j){T k ), 

J ° k=0 




where r k = 


lk/N, k = 0,l,...,N. 
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Fig. 5.3. Graphs of stresses cr n , Tn and derivatives of displacements du^/ds, du^ /ds for ctJ° = 0, cr^° = 1. 

Substituting the representations ()5.1[) , ()5.2[) into the singular integro-differential equations 
(13. lip . (|3.12p and the additional conditions (|3.10p we obtain the following system of linear 
algebraic equations: 



N 



JY 



£(<M(*i) + 9tN 2 k ( Sj )) = Re/( Si ), 5>^|( Sj ) + glNl( Sj )) = Imf( Sj ), 

k=0 k=0 

t\s)( S -l/2) k ds = 0, (5.3) 



k=0 



where sj = (2j - l)l/(2N), j = 1, 2, . . . , N, the coefficients Nl(sj), iVf (sj), N^( Sj ) and N£(sj) 
are obtained from the equations (|3.11l) . (|3 . 1 2|) . The solution of this system of linear algebraic 
equations provides the coefficients of the Taylor polynomials (|5.ip . 

Note that the system (|3.1ip . (|3.12p has additional conditions (|3.19p already built into it. 
It has been verified that the numerical solutions obtained from the system (|5.3p satisfies these 
conditions approximately. 




It has been also observed that this method provides good results even for relatively small 
values of N. For instance, fig. 15.11 shows the graphs of the function g'(s) in the case of the 
semicircular arc of a unit radius computed for N = 16, N = 20 and N — 30 and the following 
values of the parameters /i = 60, k — 2.5, 71 = 1.0. It can be seen that a reasonable convergence 
is obtained already for N = 20. Similar results have been observed in all considered cases. The 
proposed method of solution of the system (I3.1ip . (|3.12l) is computationally efficient and leads 
to numerical solution of the relatively small systems of linear algebraic equations of the order 
2N + 2, coefficients of which are computed by using numerical integration. 

Observe also that since the function g'(s) is bounded at the crack tips corresponding to 
the values of the parameter s = 0,1, it follows from the equation (13.41) that the crack has a 
sharp opening profile at the tips. This is consistent with experimental observation and is an 
improvement from the classical LEFM picture in which the crack opening profile is blunt. 

The proposed method of solution in general allows one to study an arbitrary smooth cracks. 
Some of the numerical results obtained using this method are presented below. 




Fig. 5.5. Graphs of the dependence of the maximal and minimal crack opening on the coefficient 71. 

The first studied example is concerned with a semicircular crack of a unit radius paramct- 
rically defined by the equation t(s) = e", s € [0,7r]. We consider the following two types of 
the loading applied at the infinity: a horizontal (erf 3 = 1, cr£° = 0) and a vertical stretching 
(erf 5 = 0, cr|° = 1). The physical parameters of the material are fj, = 60, k = 2.5. The depen- 
dence of the stresses and the derivatives of the displacements on the parameter s for these two 
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cases is shown on the figs. 15.21 and 15.31 correspondingly for three values of the proportionality 
constant 7! = 0.5, 71 = 1.0 and 71 = 2.0. Because of the symmetry of the crack geometry 
and the applied loading the stresses are distributed symmetrically with respect to the center 
of the crack. Observe that if the proportionality coefficient 71 = then the conditions (|3.9I) 
become (cr„ + ir n ) (so) = on the crack boundary. It can be seen from numerical results (and 
in particular figs. 15.21 and l5.3j) that (a n + «t„) ± (so) — > as 71 — > 0. 
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Fig. 5.6. Crack geometry for kq = 0.25, Krj = 0.5 and kq = 1.0. 

Here the dashed lines correspond to the values of the stresses and displacements on the 
upper surface of the crack (denoted by "+" sign), and the solid lines correspond to the values 
of the stresses and displacements on the lower surface of the crack (denoted by "— " sign). 




Fig. 5.7. Graphs of stresses cr n , r n and derivatives of displacements du-^ /ds, du 2 /ds for crj° = 1, = 1. 

It has been determined above that the derivative of displacement duf /ds and the shear 
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stress rjf may possess logarithmic singularities at the crack tips. Assume that near the tip of 
the crack corresponding to the value s = of the parameter s the functions du^ /ds and 
have the following asymptotic expansions: 

^±=A 1 \ns + 0(l), T n = A 2 lns + 0{l). 
ds 

The graphs of the dependence of the singularity coefficients A\ and A<i on the proportionality 
constant 71 are presented on the fig. 15.41 for a horizontal (of — 1, of = 0) and a vertical 
(of — 0, of — 1) stretching at infinity. The numerical results corroborate that the derivative 
du^ /ds and the stress o^ are bounded at the crack tips. 




V-i-0 

— — - y |= 2.0 



Fig. 5.8. Graphs of the dependence of the singularity coefficients A\ and A2 on the curvature kq. 



The graphs of the dependence of the maximal and minimal crack openings on the parameter 
71 are presented on the fig. 15.51 It can be seen that the minimal opening can be negative which 
suggests that the edges of the crack are pressed together shut in some areas of the crack. To 
model such behavior, one must allow for contact of the crack surfaces in the boundary conditions 
with the ends of the contact set being new parameters that must be found as part of the solution 
procedure. 

It can be seen that the dependence of the quantities in the figs. 15.41 and [5~5l on the parameter 
71 is highly nonlinear. The consideration of the graphs in the fig. 15.41 and fig. 15.51 shows that 
the extremal values of all those quantities occur for the same value of — -ff xt . The physical 
meaning and significance of the value "ff xt is immediately unclear and presents an interesting 
question to study. Since this value appears to be different for different loadings at infinity it can 
be assumed that "ff xt is dependent on the loading at infinity and not only on the mechanical 
properties of the plate material and crack geometry. 

The graphs of the dependence of the stresses o^ , and the derivatives of the displacements 
du^/ds, du^ /ds on the parameter s are presented on the fig. 15.71 for the circular arc located 
in the upper half plane between the points z = ±1 for the three values of the arc curvature: 
Ko = 0.25, Ko = 0.5 and ko = 1.0 (fig. 15. 6p . The uniform stretching (of — 1, of = 1) is 
applied at infinity. The values of other parameters are \x = 60 and 71 = 1.0. As before, the 
dashed line corresponds to the upper side of the crack (denoted by "+") and the solid line 
corresponds to the lower side of the crack (denoted by "— "). 

The graphs of the dependence of the coefficients A\ and A2 on the curvature of the circular 
arc Ko are presented on the fig. 15.81 for the following values of the parameters of = 1, of = 1 
and n = 60. Three different values 71 = 0.5, 71 = 1.0 and 71 = 2.0 of the proportionality 
coefficient 71 have been considered. 
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6. Conclusions. In this paper a new approach to the fracture modeling first introduced 
in [50], [35] has been extended to the smooth curvilinear cracks of arbitrary shape. The main 
idea of the approach is to correct the bulk material behavior in the neighborhood of the crack 
to account for effects of long-range forces from the adjacent phases. Namely, it is assumed that 
the surface of the crack is endowed with the excess properties dependent on the mean curvature 
of the crack. 

Due to the fact that the behavior of the material in the bulk is subject to the equations of 
classical linear elasticity, it is possible to apply the wide range of techniques of LEFM to solve the 
problem. In this paper the following approach has been taken. The stresses and derivatives of 
the displacements are expressed through two analytic outside of the crack functions $(z), ^(z) 
(complex potentials) with the help of Muskhelishvili formulas [TH] ■ The complex potentials can 
be further written in terms of two functions g'(s) and q(s) by using the integral representations 
of Savruk |21j . Combining these techniques with the new boundary condition which includes the 
curvature-dependent surface tension, leads to the system of two Cauchy-type singular integro- 
differential equations with respect to the functions g'(s) and q(s). This system can be further 
reduced to the system of two weakly singular Fredholm integral equations of the second type 
and, hence, possesses a unique solution except for countably many values of the parameters. 

The numerical results presented in the paper are based on approximation of the unknown 
function g'(s) by the Taylor polynomials. This leads to the solution of the relatively small 
systems of linear algebraic equations. It is observed on the examples that the proposed numerical 
approach provides a good convergence. 

It can be seen that there is a principal difference between the case of a straight crack studied 
in J23] and the case of a curvilinear crack considered here. While the stresses and the derivatives 
of the displacements do not have a classical integrable power singularity of the order 1/2 as 
in LEFM, some of them may still have a weaker logarithmic singularity. Thus, more general 
surface tension models are necessary. 
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